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We compute the two-loop heavy quark jet-function in the heavy quark limit. This is one of 
the key ingredients in next-to-next-to-leading order (NNLO) and next-to-next-to-leading-log order 
(NNLL) computations of the invariant mass distribution of top-jets at a future e'^e~ collider. The 
shape of the top invariant mass distribution is affected by large logs which we compute at NNLL 
order. Exploiting the non-abelian exponentiation theorem, a definition of the top jet-mass is given 
which is transitive and whose renormalization group evolution is determined by the cusp-anomalous 
dimension to all orders in perturbation theory. Relations of the jet-mass to the pole, MS, and IS 
masses are presented at two-loop order. 
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I. INTRODUCTION 



We are about to step into the LHC era, opening up a new energy regime for the discovery of physics beyond 
the standard model. In this era precise measurements of standard model parameters will still be important for 
disentangling new physics scenarios. A prime example of this is the effect of the mass of the top quark on precision 
electroweak constraints. The latest Tevatron analyses give mj = 170.9 ± 1.8 GeV jl), a measurement at the 1% level. 
However for a minimal standard model Higgs sector the indirect determination of the Higgs mass ttih — 76li^| GeV 
has such a strong sensitivity to the top-mass, that a 2 GeV upward shift in mt causes this mn central value to shift 
upward by 15% (with the same upward shift for the 95% CL bound mn < 182 GeV) [§. Such strong sensitivities to 
mt are also a feature of many new physics scenarios, such as supersymmetric extensions of the standard model. 

In addition to uncertainties related to the experimental analysis, the top-mass also suffers from a theoretical 
uncertainty related to the specification of the mass-scheme in which the measurement is being made. Many observables 
used in measurements of the top quark mass incorporate only lowest order theory results, which does not allow one 
to distinguish between quark mass schemes. In general quark mass schemes are connected by a perturbative series in 
the strong coupling, with relations of the form 

_schcmcA _schcmcB/i , , 2 , \ _ schemeA _schemcB i i 2 , \ /t\ 

rrif = "T-t [I + as + + . . .) , rrif = + n[as + + . . .) , (1) 

where i? is a scheme parameter (examples of schemes with relations of both of these types are discussed in the body 
of the paper) . For high precision measurements of the 6-quark mass ^ a useful class of schemes are the so-called 
threshold mass-schemes, with examples being the kinetic, IS, and shape-function schemes ^, 0, ||, 1). For the b- 
quark these schemes are optimized to avoid Aqcd sensitivities, while still maintaining a power counting in Aqcd /mi, 
to handle non-perturbative corrections in observables. The pole-mass is not used for precision analyses because of its 
infrared sensitivity, which introduces an ambiguity Sm^°^° ^ Aqcd from infrared renormalons [ p^ . For the top quark 
the infrared physics is cut off by it's width Ft = 1.43 GeV, since tops decay before they hadronize. However the top 
pole-mass still suffers from a Aqcd infrared renormalon ambiguity |ll|. 

Although we do not know the precise scheme for the top-mass measurement of Ref. jl] , we do know that it falls 
in a category of "top resonance mass schemes" , which differ from ?nf by an amount ^ F( . This follows from the 
fact that top-mass measurements rely on an underlying Breit-Wigner to incorporate the top-width, and only top 
resonance mass schemes are compatible with a Breit-Wigner line-shape [Q. For these observables using a short 
distance resonance mass-scheme avoids infrared sensitivity while maintaining a Tt/mt expansion. Examples of top 
resonance mass schemes include the jet-mass ||l^ and kinetic-mass [||. For the top-mass this scheme dependence issue 
could easily add an additional theoretical uncertainty of < 2 GeV when the measured top-mass is used in practical 
applications. To see this, lets say that the mass measurement corresponds to Tj^schcmoB^ ^^^^ we want to use a 
result for a resonance mass TTj^chomeA analysis. Since scheme B is not precisely known, the perturbative series 

introduces an additional uncertainty that we estimate to be of 0{almt) or 0{asR) from Eq. (|^), leading to the 
< 2 GeV theoretical error estimate. 

In order to obtain higher precision measurements of the top-mass one needs accurate theoretical predictions for a 
realistic experimental observable in a definite mass scheme. This has been achieved for e+e^ — > ti at threshold ||l3|, |lj, 
p^ , p^ , p^ , , where state of the art computations incorporate next-to-next-to- leading- log (NNLL) and next-to-next- 
to-leading order (NNLO) QCD corrections for the cross-section, as well as subleading lifetime effects. Theoretically 
the necessary setup is also clear for e+e" ti far from threshold, where the center-of-mass energy ^ p^ . 
Here the top quark decay products form well separated coUinear jets together with soft-radiation between the jets. 
A suitable observable is the event-shape cross-section d'^a/dM^dM^ . Here = (J^ieaPi)^ ^'^'^ ^'^i — iJ2iebPi)^ 
are hemisphere invariant masses, and the hemispheres a and b are separated by a plane perpendicular to the thrust 
axis for each event. The different physics components of d^a/dM^dM^ can be separated by a factorization theorem 
derived in Ref. ||l| 

da ,^ ,rr f Q \ f^n+^n- ( - Q^'^ ^ \o ( - 



dM^dM? 



aoHQ{Q,firn)Hra(m,—,fim,fj) [ d£^ d£ B+(st~^^,Tt,n]B^(si-^^,Tt,fj) S{£~^,£ , fj,) 
\ m /J \ m / \ m / 

+ ^ n '752^-^4'^ ■ (2) 

In Eq. (^ ctq is the tree level Born cross section, Hq and are hard-functions which encode the perturbative 
corrections at the scales Q and to, where from now on we use to for the mass of the top quark. The invariant mass 
variables St and Sf are defined as 

St M^-m^ , St Ml-m'^ 

St = — = — , St = — = — E , (3) 
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and the most sensitive region for mass measurements is the peak region where St^t ^ + QAqcd/™- Finally, B± in 
Eq. (H) are heavy-quark jet functions for the top quark/antiquark, and S is the soft function describing soft radiation 
between the jets. Our main focus in this article will be on the functions B±, which are defined in the heavy-quark limit 
rrit ^ Tt using heavy quark effective theory (HQET) |2f[ . The soft function S is universal to massless and massive 
jets and a suitable model can be found in Ref. extending earlier work in Ref. The factorization theorem 

in Eq. (j^) was derived using soft collinear effective theory (SCET) j2^, |2^, ^ and effective theory methods for 
unstable particles p^ , |2^ , A similar factorization theorem with the soft-function and different jet functions is 

known to apply for jets initiated by massless quarks |3^ . 

In this paper we carry out the first step towards NNLO and NNLL predictions for the invariant mass spectrum, 
(Pa/dMfdAI?, by computing the top quark jet function at two- loop order. We also carry out the resummation 
of large logs for this jet function at NNLL order. This translates into a resummation of all the large logs in the 
cross-section that can modify the invariant mass distribution p2| . On the conceptual side we introduce a definition 
of the top jet-mass scheme that has a well defined mass anomalous dimension at any order in perturbation theory 
(unlike definitions based on cutoff moments or peak locations). In this jet-mass scheme we prove that the quark- 
mass anomalous dimension is completely determined by the cusp anomalous dimension at any order in perturbation 
theory. As an intermediate step to demonstrating this we show that in position space the heavy quark jet function 
exponentiates. This follows from the fact that this jet function satisfies the criteria for the non-abelian exponentiation 
theorem [|6[ l^. 

Because of the simplifying nature of HQET our two-loop computation of the jet function is significantly simpler than 
a direct two-loop computation of the cross-section. In particular, as we discuss below in sections ^ and [II, even for a 
finite width and an arbitrary mass scheme the jet function computation can be reduced to the perturbative evaluation 
of a vacuum matrix element of Wilson lines. For heavy quarks two loop computations are already available for the 
partonic heavy-quark shape function |3^, ^ and heavy-quark fragmentation function |3^, ^ ^ . The hadronic 
versions of these functions that appear in observables are non-perturbative. The hadronic shape function describes 
the light-cone momentum distribution of b-quarks in a heavy B-meson while the hadronic fragmentation 

function describes the probability that a b-quark fragments to a B-meson with a particular light-cone momentum 
fraction The jet function is fundamentally different since it is defined by a matrix element evaluated between 
vacuum states, and due to the smearing from the finite top-width can be reliably computed in perturbation theory. 
We elaborate on similarities and differences below in sections || and III . 

Our outline is as follows. In section || we discuss the basic formalism for the top quark jet function, including its 
renormalization and anomalous dimension. We then give a summary of our two-loop results for the jet function and 
for the solution of its renormalization group equation, with details relegated to appendices. In section ^ we determine 
the Wilson line representation of the jet function and compare it with the shape function and fragmentation function. 
Then in section [V we work out the implications of the non-abelian exponentiation theorem for the heavy-quark jet 
function and for the partonic shape-function, including the combined implications of this theorem and the all-orders 
solution of the renormalization group equation. In section ^ we discuss possible jet-mass scheme definitions, and 
present a scheme based on the position space jet function that remains transitive to all orders in perturbation theory. 
We also give two loop relations of the jet-mass to the pole-mass, MS- mass, and fS-mass schemes. Finally, in section [v^ 
we pres ent results for the NNLO jet function with NNLL resummation, including numerical analysis. We conclude in 
section VII. 



II. HEAVY QUARK JET FUNCTION 



In this section we describe the basic properties of the heavy-quark jet functions B±. Up to a change of variable 
-B_i_ for the top quark and B_ for the antitop quark are identical by charge conjugation, so we will only refer to the 
computation of B+. To simplify the notation we also drop the subscript, so that B = 5+. These subscripts ± are 
restored when we simultaneously consider the top and antitop system in the final factorization theorem. We start by 
reviewing definitions and results for the HQET jet function from Refs. |T^, |T9j. B is given by the imaginary part of 
a forward scattering matrix element, 

B{s,6m,rt, fJ.)^lin[B{s,Sm,rt, fi)] , (4) 
where B are vacuum matrix elements of a time-ordered product of fields and Wilson lines 

B{2v+-rJm,Tt,^i)^-^- /d^a; e''-" ( |T{/i, jO)W„(0)M^^(x)/i„^ (a;)}| O) . (5) 

Here is the velocity of the heavy top quark, and we introduce null-vectors n'^ and n'' so that we can decompose 
momenta as p'^ = + n^n-pjl +v\- The vectors satisfy v\ = \ and = = 0. The dot-products of these 
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vectors encode the boost of the top quarks relative to the center-of-mass frame of the e+e collision, n ■ v+ = m/Q, 
and n ■ v+ — Q /m. In Eq. ^ the Wilson lines are 



W^(a;) = P exp J dsn ■ An{ns + x)^ , Wn{x) — P exp ^9 J dsn ■ An{ns + x)^ 



(6) 



These Wilson lines make B gauge-invariant and encode the residual interactions from the antitop jet. Both the HQET 
fields hy^ and the gluon fields in Wn (which we call A^) are only sensitive to fluctuations with <C m^. In the rest- 
frame of the top quark these are soft-fluctuations, while in the e~^e~ center-of-mass frame they are "ultra-collinear" 



along the direction of the energetic top quark. The gluon fields A!^ have zero-bin subtractions |45 for the region of 
the soft function S in Eq. as explained in Appendix B of Ref. . 
The HQET fields h^^ have the leading order Lagrangian 

Ch = ■ D - Sni+ ^Tt^ hy^ . (7) 

Here Ft is the top quark total width, obtained from matching the top-decay amplitudes in the standard model (or a 
new physics model) onto HQET at leading order in the electroweak interactions, and at any order in a^. This gives 
the correct description of finite lifetime effects for cross-section in Eq. (H) to 0{m^ /Q^ ,r /m) in the power counting 
for separation of the jets from the decay products ]T^ . The residual mass term 5m in Eq. (|^) fixes the definition of 
the top mass m for the HQET computations where 

Sm^mpoie-m. (8) 

For predictions in the peak region consistency with the power counting requires Sm ^ Ft ~ st ~ fl^ , a condition 
which is true of the jet-mass scheme that we discuss below in section ^ 

From the definitions in Eqs. (|^^ and the Lagrangian in Eq. (|^) one can deduce a series of properties of the 
jet function. As a first, instead of computing _B(s, (5m, Ft, /i) and B{s, 6m,Tt, n), one can consider computing these 
functions for a (fictitious) top quark having zero width. Furthermore, due to Eq. (^ the s and Sm dependence occurs 
in the combination (s — 2dm), so it is useful to also have a notation for computations done with a zero residual mass 
term in the Lagrangian. Thus we define 

B{s, Sm, fj.) = B{s, Sm, 0, ^) , B{s, Sm, /i) = B{s, Sm, 0, fi) , 

Bis,fi) = Bis,0,0,fi), Bis, fi) = B{s, 0,0,ii). (9) 

These jet functions and vacuum matrix elements are related by 

B{s,Sm, fj,) — lm[B{s,Sm, n)^ , B{s, fi) — lm[B{s, fi)] , (10) 

and B(s,fj,) has support for s > 0. The form of the Lagrangian in Eq. (^ implies that having calculated B{s,fi) we 
can include the width and Sm terms by simple shifts, 

B{s, Sm, Ft, /i) = B{s + iVt, Sm, /i) = B{s — 2Sm + iVt, /i) . (11) 

As discussed in Ref. |l^ the stable and unstable HQET jet functions can also be related with a dispersion relation, 

Tt 



ds' B{s -s', Sm, fi) ,.,n\ ^ / ds' B{s ~ s' - 2Sm, fi) 



12) 



The width of the top quark acts as an infrared cutoff through this smearing with the Breit-Wigner. Finally we remark 
that the /i-dependence indicated by the last argument of B{s, Sm, Ft, /i) and B{s, Sm, Ft, fj.) is independent of Ft and 
Sm. Additional scale dependence may be induced by the choice of mass-scheme, ie. by a parameter Sm = Sm^jj). 
When we consider B{st,Sm, Ft, fJ,) as a function of Mt this additional /i-dependence from Sm cancels against that in 
the mass m{fi) in Eq. (0). This cancellation occurs at leading order in the HQET power counting. 
We will also find it useful to consider the Fourier transformed jet functions 



ds e-'y^ B{s,Sm,Tt,fJ.) , B{y,Sm,fi)= ds e'^^ ^ B{s,Sm, n) , (13) 

-OO J —OQ 

where y = y — iO to ensure convergence as s — s- oo. In Fourier space the connection between the jet functions computed 
with zero and non-zero width and residual mass terms becomes particularly simple, 

Biy, Sm, Ft, Ai) - B{y, Sm, fi) e"!^! = B{y, fi) e-"^^*'™ e"!^! . (14) 



This formula is quite interesting, since as we discuss in section IV below, the result for B{y,ii) also exponentiates to 
all orders in perturbation theory. 
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A. Renormalization and Anomalous Dimension 

We use dimensional regularization with d = 4 — 2e and the MS scheme to renormalize the jet function. The 
renormahzation properties of B(s, fx) and i3(s, n) are the same, so in the following we work with i3(s, yu) for simplicity. 
The divergences of loop calculations are removed with Z-factors, so one can pass from bare to renormalized matrix 
elements by 

S(J,^)= J ds' Zg\s- s',^i) B^''"'is'). (15) 

This equation can be thought of as the generalization of a Z matrix which renormalizes a set of operators indexed by 
s, to the case where s is continuous [M. Here Zb and its inverse satisfy 



ds Zg^{s" — s ,^)Zb{s — s,fi) = S{s" — s) . 



From the /i independence of one obtains the renormalization group equation 



d f 



where the anomalous dimension is 



/d f d 

ds' Zg^{s - s',fj.)fi—ZBis',^J.) = J ds' Zb{s - s' , ^i)^i—Z]^^ {s' , ^i) 



(16) 



(17) 



(18) 



Since ^b{s,^) is real we can also simply take the imaginary part of Eq. (|T^) to obtain the renormalization group 
equation for B{s, /i). In the MS scheme Zb and Z^^ have the e dependence 



oo 

ZB(J,Ai) = '5(s)+5]-ZW(J,M), 



OO 

Z5i(J,M) = '5(5) + ^-ZW(J,M), 



(19) 



k=l 



k=l 



where ZC'^) and Z'*^) are e independent. Eq. (H) implies that Z^^) = -Z^'^\ and Z^^' = -Z^^) - Y.]J^ Z^-?) ® Z^^' 
for k > 2. Demanding that 7_b(s, /i) is finite as e — s- and using the /3-function equation 



^l—as{^J.)^-2eas{^i) + f3[as] , 
d^ 



(20) 



gives the standard dimensional regularization result that the anomalous dimension is determined by the residue of 
the 1/e term at any order in perturbation theory, 



7s(s,m) = -2a, — Z«(s,//) = 2a,— Z«(s,//). 



(21) 



We find that the higher 1/e poles lead to the consistency equations [i > 1] 

2a.^4'^^'(5,A^)-(M|-+/3^)4'^(^,M) + E Jdrs' zPis-s',^) 

k—1 



2a. 



das 



+ M 



Ofi oas ) 



(22) 



where for convenience we let Z^^(s^[i) — b(s). The result in Eq. ( [2^ ) agrees with the form of the counterterm 
consistency condition derived in Ref. |Q for the heavy-quark shape function. 

At any order in perturbation theory the anomalous dimension in Eq. ( pT| ) has the form 



7B(s,M)=-2rK 



^iQ{s) 



7[as]5(s) 



(23) 



where our definition of this plus-function is given below in Eq. ([10[). Here r'^[as] and 7[q;s] have an infinite power 
series expansions in a, that starts at linear order. r'^[as] is the cusp- anomalous dimension p7| , 49, 5^, while 7 [a,] 
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is the part of the anomalous dimension that is unrelated to the cusp. In position space the renormalization group 
equation and anomalous dimension are simpler, 



lB{y,f^) = 2r'^[as] In(ie^^y^) +7^ 



(24) 



The form of the anomalous dimensions given in Eqs. ( |23|) and (|24| ) is guaranteed to all orders in perturbation theory 
by a theorem regarding the renormalization of Wilson-line operators with cusps proven in Ref . Eol . pO| , which ensures 
it can not have dependence on the position space variable other than the ln(j//i). To solve Eq. (p4|) one first writes 

ln(ie'''^y/Lt) = In^ie^'^yno) + ln(/i//io), then rewrites ln(^/^o) = da' /f3[a'], and finally integrates with a change 

of variables dln/x = da/P[a]. This gives a solution that connects the result at the scale /Xq to that at the scale /Lt, 



(25) 



where the two evolution functions are 

da 



w(Ai,/io) 



/3[c 



■T^[a]. 



da 



7N + 2 



"^^'^^ da 



da' 



(26) 



Taking the Fourier transform of Eq. (^5|) then gives the solution to the momentum space renormalization group 
evolution (RGE) equation 



B(s,m) 



ds' Usis - s',iJ.,iJ,o) B(s',no) , Ub{s - s',iJ,no) 



(s - s')i+" 



(27) 



where K = K^h^hq) and lo = ui{ijl, Hq). All results in this subsection are valid to all orders in the as expansion, 
and can thus be used to sum logs in B at leading log (LL), next-to-leading log (NLL), NNLL, and beyond. To our 



knowledge, the results in Eq. ( |26D and 
Ref. and then to all-orders in Ref. 



were first derived for the i?-meson shape function, first at one-loop in 



B. NNLO Result for S(s,^) 



To obtain results at NNLO we consider the as expansion of quantities defined in subsection 
renormalized jet functions can be written as 



[I A. The bare and 



J=0 



6°"'-'= (J) 



B 



3=0 



We also expand the anomalous dimensions and /3-function as 

-a,(^)iJ+i 1 ras(M)i^'+i 



r[as]=Y.T^[- 



An 



3=0 3=0 
where up to three-loop order [|3[ |^ ^ ^ | 

3~' 



47r 



(5[as] = -2a,(^)^/3„ 



n=0 



as{^i) 



Aw 



Tl+l 



/?2 



2857C' 
54" 



18 



54 



9 



3 

•9C 



(28) 



(29) 



(30) 



To incorporate the (5m term from the Lagrangian in Eq. (|^) we evaluate ;B^''"(s — 25m) and then expand in as(/i) 
with 



(31) 
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This is simpler than treating Sm as a Feynman rule insertion, and equivalent. The bare and renormalized couplings 
are related by 



a 



bare 



t = exp(7B)/(47r) , 



(32) 



where Zg is the Z- factor for the strong coupling and the iota dependence, l'^, ensures we are in the MS scheme rather 
than the MS scheme. To determine the renormalized jet function we expand the counterterms as 



on , , 



1 l-si,)V^i^^^,^^^^ 



k=l j=l 



(33) 



1 



Using this notation, converting a^^''° to as(^) with Eq. (|3^), and then equating powers of as(M) Eq. ( p^ ) these 
expansions determine the renormalized Bj(s,dm, The tree, one-loop, and two-loop coefficients are respectively, 



6o(s,5m,Ai)=6o^--(s), 



(34) 



dBo{s,fj,) 
ds 



B2{s,5m,^i) = ^'^x'''B^"{s) + 2zgli'^?'B^''{s) + / ds' Zi{s^s' , ^i) L^fi^'Br^is) + / ds' ^2(5-5', /i) Sj^""(s') 



dBois,fi) , , 
2 0TO2 77 1- 2 (omi) 



2 d'^Bo{s,fi) 
ds ' ^ c?^s ds 



— 2 (5toi 



dBi{s,fi) 



where we used a subscript notation for the [as(/i)/7r]-' expansion coefficients as in Eq. (^). The one and two-loop 
Z-factors have terms 



Zi 



1 



where the coefficients zj'^'' are defined so that the Bj{s, Sm, ^) are finite as e ^ 0. 

The results for Bq and Bi were obtained in Ref . [^2|, |T9j . In an arbitrary mass-scheme we have 

mBo{s, 5m, fi) ^ L" , m6i(s, Sm, fi) = Cf[l^ + + (l + ^) L°} - 2 Jmi(L°)' , 

where the prime denotes a derivative with respect to s, and for convenience we have defined 



1 



^(2) 



z. 



1 



1 



^(2) 



1 



,2 ^2 



^(3) 



1 



,3^2 



(35) 



L'--^ . ^ ln^-(- 



7r(— s — iO) 



s - iO 



(36) 



(37) 



The corresponding two-loop result, B2, is one of the main results of this paper and involves the Feynman diagrams 
shown in Fig. |l|. Details of the computation of B2 using Eq. (|34|) i n Feynman gauge are given in Appendix To 
summarize, we use the computation of Broadhurst and Grozin |59[^0| for the divergent and finite terms of the two- loop 
heavy quark propagator (the first graph in Fig. |l|), and compute the remaining Feynman diagrams directly. We treat 
the quarks other than top as massless, with n/ such flavors, and thus do not include effects due to the 6-quark mass 

~(k) 

in vacuum polarization diagrams. We have also confirmed that the resulting Z^ satisfy the counterterm consistency 
conditions in Eq. (^). The final result for the renormalized two- loop matrix element is 



771^2(5, Sm, n) — C 



1 



2 

C'fCa 



1 TT 



137r^ 
5 



1 H 4C3 i + - H 

24 ^'V V2 24 

11 Stt^ IOtt" 



537ij 
"640" 



-2C3U'' 



V3 12/ V18 12 4 / V 54 48 



960 



5C3 
8 



2.0 /47 TT^x 1 /281 237r2 17(3 
6" '3" ■ \^^\2) ^V216^T92 48" 



2Sm2{L")' + 2{SmiY{Ly - 2Smi Cp 



5^ 
'24 



(38) 
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FIG. 1: Graphs for Gluons from an (g) are from the Wilson hnes, the hatched blob is the two-loop vacuum polarization 

of the heavy quark, and the blobs with diagonal lines include all one- loop vacuum polarization graphs for the gluon. Numbering 
the graphs from 1 to 16 from left-to-right and top-to-bottom, we note that graphs 2, 4, 5, 7, 9, 10, 11, 14, and 16 implicitly 
also stand for their left-right symmetric counterpart. 



One can pass from the funetion 6(s, dm, n) to the distribution B{s, dm, ji) using Eq. ( [To| ) and the formulas in Eq. 
This gives 



m B{s, 5m, m) = 5{s) + ^^I^M Hc'-C'^+il- '- )5[s) 



2£^ - 



11^2 1 

3- — — 



2 a sin) 



(5mi(/i) S'{s) 




12 

5 7r2 5C3 

\ V — 

18 12 4 



24 



1 Svr^ IStt^ 



2^3^V 36^12^^1216 576 48 



24 5760 
54 144 2880 8 / ^ _ 

2 



2C3M(s) 



(5m2 (5'(s) - {8mxY (5"(s) + (5toi Ct 



2{C')' -{C 



0\/ 



(39) 



where for the log plus-functions we use the notation 



-I + 



lim — 



In 2^ 



In'^+iC 



(40) 



x—sj 11 



The results in Eqs. (|38|) and (^9|) are presented in an arbitrary mass scheme, which is specified by the choice for the 
coefficients 8m\ and 5mi. An appropriate mass-scheme for top-jet cross sections is described below in section 0. 
In order to obtain the distribution B{s,Sm,rt, n) with Ft one can input Eq. ( p9|) i nto the integral with the 
Breit-Wigner in Eq. (^2|). However the simpler method, which we use below in sectioii"|v^ is to shift s ^ s + iTt in 
Eqs. ( |36| , p8| ) and then take the imaginary part as in Eq. (pi]). 

From the renormalization constants z[^^ and ^2^' given in Appendix ^ we also obtain the anomalous dimension 
terms in Eq. ( p9| ) . The cusp anomalous dimension is known up to three-loop order Fq 2 JSOl , |6l| , and we have 
confirmed that our two- loop analysis reproduces the expected result for the two-loop cusp coefficient, F^. For B the 
one-loop anomalous dimension 70 has been calculated in Ref . |19] . The two- loop anomalous dimension 71 is obtained 
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from our calculation of Z^^-* in Eq. (All).^ We list here all the pieces needed for our analysis, 



4C 



/268 47r2\ „ „ 40 



V 9 

490 5367r2 447r4 
~3 ^ 



CfCa — Tr^pnf, 



27 



45 



-836 112C3 



70 = 4Cf 



71 



1396 237r^ 



27 



20(3 



27 

CfC'a - 



2^ 
9 



232 
"27 



CprifCA 
Cprif . 



32C3 



110 



WCpnj 
27 



(41) 



To resum the large logs in the jet function to NNLL order we need to use these coefficients in the NNLL results for 
the evolution functions uj{fi,fio) and K{fi,fio)- To NNLL order solving Eq. (p^ ) gives 




/3i\ (1-r+lnr) /3i 



327r2 



in 



In^ r 



(/3o7i-/3i7o) 



Po^o Pq 



P2 PI 



^ (l-r+rlnr)+ ^ - ^ (r-l) Inr - 



Po Pt 



PiVl P2 
Pq^o Po 



167r2 

/3p(l-r)2 
- 



(42) 



(r-1) 



P? 



where r — as{iS)/as{fJ-o). Eq. ( ^ ) determines the evolution functions in terms of coefficients of the anomalous 
dimensions and /3-function. It agrees with the NNLL result given in the appendix of Ref . , which was used to sum 
large logs in the _B-meson shape function for B Xg'j- 



We postpone presenting our final resummed NNLL result for the jet function until section VI 



III. WILSON LOOP REPRESENTATIONS, AND COMPARISON WITH THE HEAVY QUARK SHAPE 

FUNCTION AND FRAGMENTATION FUNCTION 

It is well known that the leading order coupling of gluons to heavy-quark fields hy in HQET can be represented by 
Wilson lines along the path of the heavy-quark |6^. We define 

Wy{x) =Pexp(^- igj ds v-A{vs+x)J , W^{x) = P exp (^ig J ds v-A{vs + x)J . (43) 

To see how the HQET action reduces to a Wilson line one can make a field redefinition, hy — Wyh^\ from which we 
find that hy ' is a free field with Lagrangian Ch = hi, iv • , see ||26| . Thus, the vacuum matrix element for the 
heavy-quark jet function in Eq. (^) can be written as a matrix element of Wilson lines 

S(2t;-r,^) = — ^ fd''xe''---{0\Th^°\0)W^{0)W„iO)W,lix)Wyix)h(^Hx)\0) 



2'ENcra 
i 

2TTNrm 



dx° e™-"° e{x°) { \tTTW^^{0)Wn{0)Wlix°)Wyix°)\0) 

e™-"° 6'(x°) (0|trTW^(a;",0)VK„(0,oo,x")|0) , (44) 



where 2v-r — s and we use the shorthand x'^ — v-x, the trace tr is over color indices. Here Wy{x'~',0) — Wy{x'~')W^{0) 
is the straight Wilson fine from to while W„(0, oo, a;") = W^„(0)W^(x°) has a path from x^ to oo to that 
uses two light-like line-segments. To obtain the second line of Eq. (^) we used the heavy-quark propagator, 

{0\Th^^''{0)h^y^\x)\0) = -26''''S^ix)0ix°) where a and b are color indices. In Fig. Ha we give a graphical repre- 
sentation for the Wilson line definition in the last line of Eq. (H) . The arrows denote the time-ordering. 



It turns nut that the piece of 71 proportional to Cpn.f is the analog of a contribution in the analysis of a scalar field theory made in 
Ref. pa, pu]. Suitably translated to the QCD case their computation agrees with the Cpnf term in our 71. The non-abolian CpCa 
term of 71 is original to our work. 
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FIG. 2: Representations of the Wilson line matrix elements for the heavy quark jet function, The Wilson lines include Wv, 
(double lines), and W„, WX (single lines). In a) we display the result in Eq. (^) that gives B{s, fi). In b) we display the result 
in Eq. (^ ) that gives B{s,fi). 



We can also write the jet function B{s, /i) as a matrix element of Wilson lines. To derive this result we note that 
B{2v■r,^i) = —J— /d^xe*'--" V Tr {0\TW^{x)K{x)\X){X\Th40)Wn{0)\0) 



1 



d'^x e^"--" Tr (0| [TW^{x)K{x)] [Th,{0)WniO)] |0) 



= Jd''xe^^--TT{0\[TW^,{x)W4x)h<;^\x)] [Th'^\0)W^{0)WniO)]\0), (45) 

where T is time-ordering, T is anti-time-ordering, and the trace Tr is over spin and color indices. Next we use 

(0|ft.[,^^°(x)ft.[,°^''(0)|0) = 25"^(5^(x), where here there is no time-ordering and hence no 0{x°), and the spin indices i are 
contracted. Thus 

B{S,fi) = ^ fdye^'yBiy,ti), B{y,^i) = -^{o\tT[TW^{2y)W42y)][TW^iO)WniO)]\0). (46) 



Here we took x^ — 2y in order to agree with the notation for the position space jet function B{y, fi) above in Eq. (|lJ 
In Fig. ^ we give a graphical representation for the Wilson line matrix element for B{s,ii) in Eq. (|4^), where the 
arrows denote the time-ordering. Comparing to the Wilson loop for B{s,fi) in Fig. ||a we note the importance of the 
oo-points to determine which fields are time-ordered and which are antitime-ordered. 

It is instructive to compare the Wilson line definition of the heavy quark jet function with the corresponding defini- 
tions for the heavy quark shape function that appears in i?-decays |4^, 43 , and with the heavy quark fragmentation 
function [Q. Using a variable ^""^ > 0, the B-meson shape function is given by 

f''-ii+,f^) - ^ fdx~ e-''(^^-^)^"/2 J2 {B.\TK{0)Wn{0)\x){x\TWl{x)K{i)\B,) 



X 



= ^ (dx- e-(^+-A)--/2 (5^1 [tU^\0)WI{0)W^{0)\ [TWl{x)W,{x)h'^^\x)\ 



87r 

^ ^ jdx- e-^(^+-^)--/2 (B„|/i(o)(0)VKj(0)I¥„(0,i)ll^,(£)/i(")(i)|B„) , (47) 

where we use the shorthand x~ = n - x and x^ = x~n^ /2, and set n ■ v = 1. To obtain the second line we made the 
same field redefinition as above in Eq. (^), but now on both the heavy-quark field and on the interpolating field for 
the i?- meson [^3[ Q . Combining the lines from these sources yielded the Wilson lines 

Wy{x) ^ P exp iig i ds v- A{vs+x)\ , W^{x) exp ( - ig dsvA{vs + x)j. (48) 

To obtain the third line of Eq. (g^) we noted that TW^ = W^, TWn = W„, TW^ = Wy, and TW^ = W'J, and that 
the gluons in the Wilson lines which sit next to each other, WX{x)Wv{x) and W^J(0)H/Vi(0), are already time-ordered 
and anti-time-ordered respectively. 

For the B-meson fragmentation function in HQET with variable > 0, the field redefinition gives 

D''^/\e+,fi) = fdx- e»(^++^)-"/2 ^ (0|rW^,U^)/i.(^)|^.X)(B,X|T/;,,(0)W^„(0)|0) 

" 1^ Z'^^" E (0|rW^,Ui)W^.(5)4°H5;)|^.X)(B,X|T;i(o)(0)W^J(o)M/„(0)|o) . (49) 
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Here Wn and are defined as in Eq. (|8|) but with v ^ n. The shape and fragmentation function results in Eq. ( |47| ) 
and ( ^9| ) are similar to the heavy-quark jet function in that all three are defined by matrix elements with heavy-quark 
fields and Wilson lines. They differ because they are non-perturbative distributions involving a B-meson state in 
contrast to the perturbatively computable jet function. The shape and fragmentation functions also have a light-cone 
separation rather than the time-like separation that we have for the jet function. 

In certain contexts it is also useful to consider the partonic shape function Z**" and the partonic fragmentation 
function D^"/'' where the _Bt,-meson state is replace by a 61,-quark with residual momentum A:^, where A = 0. In this 
case we can perform the contraction hll^\x)\by) = e~*'°'^u„|0)/\/]Vc and write 



1 



dx e 



-i(£++k+)x- /2 



(0|trWj(0)W„(0,£)iy„(£)|0) 



(50) 



Attn, 
1 

1 



-Nr .1 \ 

dx- e*(^+-'=+)-V2 (^o\trWl{0)Wn{0,x)W,{i)\0) 



where we used Wn{0)W^{x) = Wn{0)W^{x) = Wn{0,x). Eq. ( |50| ) states that the partonic shape function and 
fragmentation function are identical, but with complementary ranges of support. This was observed in Ref. p9| for 
logs at NNLL accuracy, and was derived to all orders in perturbation theory in Ref. |^ as we outlined above. Thus, 
the partonic shape function in position space, {y, ^) is also given by a vacuum matrix of Wilson lines. It differs 
from B{y,jj,) in Eq. (^ ) both due to the light-like rather than time-like separation y — x- /2, and due to the path. 



IV. NON-ABELIAN EXPONENTIATION 

In the previous section in Eq. ( ji^ ) we showed that the position space heavy-quark jet function B{y, /i) is determined 
by a vacuum matrix element of Wilson lines. Due to the non-abelian exponentiation theorem [ p6| for matrix 
elements of Wilson lines with symmetric restrictions on the phase space of real gluons, /x) exponentiates. This is 
also true of the partonic heavy-quark shape function in position space. 



Taking the Fourier transform of the two-loop jet function result in Eq. ( |39| ) using Eq. (Bl) we obtain 



mB{y,fj,) = 1 



TT 

24 



CfPo 



1 ~, 
6^ 



36 48 



57r2 
576 




11 

54 



281 
216 



(51) 



where L*"' = (L)'' and 



L = In {ie^'^y fij 



(52) 



It is evident in Eq. (^T|) that the two-loop, Cpa^ term satisfies the exponentiation theorem, being determined by 
one- half the square of the one- loop Cpots term. Thus we can write 



mB(y,^) = exp|^^:^^lMfp 



TT 

24 



o^MCfCa 



V3 12/ 



' 18 12 • ' 



5C3 



2 - 



36 48 



177r4 
2880 



576 


+ 7 


7^2 


11 


144 


54 



(53) 
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The non-abelian exponentiation theorem guarantees that corrections to this result are 0{a^) in the exponent, and 
that these corrections will vanish if we take the abelian limit Ca — *■ and — > 0. 

In the abelian limit with zero /3-function, the exponentiation theorem implies that ln[mi3(y,/i)] is one-loop exact. 
Thus taking Ca ~ rif — 0, and a charge Cp we have to all orders in perturbation theory 



abelian 



exp 



in 



- 7o£ + To 



(54) 



where the constants are Fg = 70 = 4Cf and Tq = 4Cf(1 + 7r^/24). The exact result in Eq. ( p^ provides a simple way 
of testing the properties of different possible jet-mass definitions at higher orders in perturbation theory, as discussed 
in section ^ . 

We can a so consider the implications of the non-abelian exponentiation theorem for the solution of the renormal- 
ization group equation (^). Following Ref. we first use the evolution kernel K[ij,,^q) to solve for B{y,iJ,) by 
taking /hq = Hy 



-IE 



jy. This makes all the logs in B(%j^ \iy) vanish since L{fjLy) = \n(ie'^'^yiiy) — 0. Thus 
mB{y,fi) = e'^^f^'f^y^ mB{y,fiy) = eK{^.,^.y)+T^o.,{^y)] 



(55) 



Here the boundary condition for the RGE, denoted ■ml3{y, fiy), is just a perturbative series in Usiny). Due to 
the non-abelian exponentiation theorem this series must exponentiate to give exp(T[as(/iy)]), and the coefficients 
in the perturbative series for T[as] have color factors that satisfy the exponentiation theorem constraints. It is a 
straightforward exercise to verify that expanding the result for _fC(/i,/io) in Eq. (^) to 0{a1) gives a result from 
Eq. ( p5| ) that is consistent with Eq. (p3|). 

The Fourier transformed partonic b-quark shape function is also given by a vacuum matrix element of Wilson lines 
via Eq. (|5^). Thus, it too satisfies the criteria of the non-abelian exponentiation theorem p8[ . Taking the Fourier 
transform of the two-loop computation of f^^'{t^,^) in Ref. |Q we have verified that the C^a^ terms satisfy the 
non-abelian exponentiation theorem. This calculation gives 



{y-, m) = exp ■ 



+ 



'24 
12/ 



a2(^)CF/3o 



11 
18 



TT 

12 



6 6 



9C3 

4 



9C3 



(l-37r^ 
36 

IOTtt^ 



-L- 



2880 



IIC3 



48 
13^2 



III 
192 



29 
108 



1 
216 



(56) 



where L is defined in Eq. (|52|), but now y is the conjugate variable to £^ , y 



/2. Corrections to this result are 



again ©(af) in the exponent, and vanish when Ca — nf — 0. Comparing Eq. ( p3[ ) and ( |56| ) we explicitly observe the 
difference between the heavy quark jet function and the partonic shape function. Up to a sign the highest powers of 
L agree at each order in as , because of the relation between their cusp anomalous dimension terms. The subleading 
logs and constant terms differ. 



V. A TRANSITIVE JET-MASS SCHEME 



The last remaining ingredient needed for the NNLO and NNLL computations of the heavy-quark jet function is the 
specification of the mass scheme counterterm Sm at two-loop order. Since the jet function will be used to describe 
momenta s ^ F, where F is the width of the physical invariant mass distribution, we must have 5m ~ F or smaller to 
not upset the power counting in the HQET Lagrangian, Eq. (^. In the MS scheme SrfT ~ m{as + + . . .), and since 
mas 3> F this scheme does not satisfy the power counting criteria. In the pole-mass scheme Sm = to all orders, 
however this scheme has instabilities related to its infrared sensitivity. In particular the pole-mass has an infrared 
renormalon that leads to an asymptotic ambiguity Sm^"^'^ ^ Aqcd7 and hence is not a useful scheme for precision 
computations. Schemes that satisfy 5m ~ F and do not suffer from infrared renormalons were called top "jet-mass" 
schemes in Ref. |T^. We refer to them more generally as "top resonance mass-schemes" here and reserve the name 
jet-mass for a specific example of this type of scheme. These mass-schemes are suitable for use in the factorization 
theorem for the top-invariant mass distribution in Eq. (3) and related observables. We start by defining a jet-mass 



mass schemes in section VB 



scheme with nice renor maliz ation properties in section [V A|, and then relate this jet-mass to the pole, MS, and IS 
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Potential Jet-Mass Definitions and Anomalous Dimensions 



by 



In this section we explore three resonance mass-schemes for m. With the notation for dm in Eq. (H) they are defined 



a) 
b) 



ds 



B(s,^mP<=^^^t,^i) 



= 0, 



(57) 



f ds s B(s,(5m'"°'°,^) = 0, 

J — OO 



2B{y,ij) dy 



y = -ie-1E /R 



^ R d 
2 d\n{iy) 



iye1E=l/R 



We refer to a), b), c) as the peak-mass, moment- mass, and position-mass respectively. The peak- mass definition uses 
the jet function with a non-zero width and satisfies the 6m ~ Tt power counting criteria |p^ . In b) and c) the schemes 
depend on a parameter R, and we must take i? ~ Ft in order to satisfy the power counting criteria. Different choices 
for R specify different schemes, and are analogous to the difference between the MS and MS mass-schemes. All three 
schemes in Eq. ( p7[ ) are free from leading renormalon ambiguities |65 . In the following we will argue that only the 
definition in c) is a reasonable scheme for higher order computations. Thus we will only use the name jet-mass for 
this position-scheme mass definition. 

The definitions in Eq. ( |57| ) are all perturbative mass-schemes which stabilize the peak position of the jet function 
B{s, Sm,Tt, ^i). In scheme a) the peak position is fixed to all orders in perturbation theory by definition. In scheme 
b) we instead fix the first moment, which provides a more local observable that is still sensitive to the peak location. 
However, scheme b) still has non-locality induced by the cutoff R on the momentum space moment. A finite R 
is necessary due to ultraviolet divergences that occur for i? — > oo. This type of moment divergence is a general 
property of functions that have a cusp anomalous dimension (see for example Refs. [|66t ^7\). If it was not for the 
UV divergences then the schemes b) and c) would be equivalent in the limit R ^ oo. In the situation at hand, c) 
provides an independent mass scheme definition. A jet-mass definition from c) is explicitly local since it just involves 
the position space jet function at a particular position y. 

An additional criteria for a reasonable jet-mass scheme is to have a renormalization group evolution that is transitive, 
as discussed in Ref. [19 . Transitivity is a well-known feature of the MS mass, and implies that we will obtain the same 
result if we evolve directly from /iq — > /i2, or if we first evolve from /io — > /ii and then from /^i — > //2- Transitivity is 
guaranteed by any mass-scheme with a consistent anomalous dimension and renormalization group equation. Since 
in HQET the scale independent m^°^° = m{fj,) + Sm{fj,), the general form for the RGE equation for the mass is 



7m = -^J■-l-^'m{^J-) , 



(58) 



where i? is a mass dimension-1 scheme parameter. Transitivity of m(/i) is guaranteed by this anomalous dimension 
equation, as long as -fm is proportional to [m{ij.)]'^ R^~^ for some k (and thus, for example, is not a sum of two types 
of terms with different powers of fc). In the MS scheme k = 1 and the anomalous dimension is proportional to m{fi), 
while in all three schemes in Eq. (57) we have fc = 0. However, it turns out that the peak-scheme and moment-scheme 
do not have consistent anomalous dimension equations of the form in Eq. (|5^), because there 7,„'s depend on explicit 
powers In^ (fi/Tt) and In-* (/i/i?) with higher and higher powers of j > 1 occurring for higher orders in a^. These logs 
render the moment scheme anomalous dimension equation inconsistent at NLO order, and the peak scheme does not 
even have an anomalous dimension equation of the form in (|58| ) at LO order. 

In order to illustrate the difference between the three schemes in Eq. (^^ we first consider the simplified case of 
the jet function in the abelian limit, Ca and n/ 0. The all-order result for B{y,fj,) is given in Eq. ([m|), 
and can be directly used to determine Sm in the position-mass scheme. The derivative of the exponential gives back 
an exponential which cancels against the 1/B{y,fi) in Smj. Thus the abelian result in the position-mass scheme is 
one-loop exact. 



R 2 



(59) 



Since for the abelian limit dus/d^ — 0, the abelian anomalous dimension computed from Eq. ( [59| ) is (7;^)'*^''''™ = 
—Re^'^CFOis/i^ to all orders. Thus this position-scheme anomalous dimension has the desired form in Eq. (p8|). To 
compute results for the peak and moment mass-schemes we need the abelian jet function in momentum space, B{s, /i). 
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Tree, one- loop, and two-loop terms are given by the abelian terms in Eq. To determine three-loop and higher 

order results we can simply expand Eq. (^ ) in and take the Fourier transform. For the three-loop term in the 
abelian jet function this gives 



TOi?3(s,/i) 
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2 



abcl: 



lan 



£5 - 

2 



IStt^ 257r4 
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(60) 



^^ + ^-12C5)<S(s) 



we find that up to three-loop order 
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(61) 



At one-loop order the Sm factors in the three schemes a), b), c) are quite similar. ^ However the three schemes are 
quite different at two-loop order. Computing -fm from these counterterms we see that a ln(/i/Ft) appears in 7m°abciian 
at two-loop order, and that a ln(^/i?) appears in 7m°Jboiian three-loop order. At higher orders in a^, higher and 
higher powers of these logarithms, ln-'(^/i?), appear in 7„i in the peak and moment schemes. In particular we see 
from Eq. (|6^) that at three-loops for the peak-scheme there is a Cp ln'^(/i/Ft) term that generates a \n^{^/rt) in the 
computation of 7^- For the moment scheme we have extended the abelian computation to four- loops, and find a term 



{Sr 



\ 41oop 



abelian/ 



/32 47r4\ 

vy ~ I5"/ 



In 



3M 

R 



(62) 



This gives a \n^{fi/R) in the moment scheme 7„i at four-loop order. ^ The absence of a^' In'^ (/i/i?) terms in ^TiJJ^^J^j^ is 
a reflection of the fact that the moment- mass has a consistent anomalous dimension at LO. Neither the peak-scheme 
nor the moment-scheme have consistent anomalous dimension equations in general. This inconsistency arises because 
of the non-locality inherent in their definitions. This is problematic because we would like to be able to evolve our 
mass as a function of /i, for instance to run it up to large mass scales and connect it to the MS scheme at a scale 
/i = mf. On the other hand the position-scheme is entirely local, and so far we have demonstrated that in the abelian 
limit it has a consistent mass anomalous dimension. We extend this proof to the non-abelian case below. 

In table || we present non-abelian results for the two- loop computation of 6m for all three schemes in Eq. ([5^). 
The position-scheme jet- mass is no longer one- loop exact, and has corrections at each order in as{fi). Since now the 
coupling as{n) evolves, higher powers of In-* (/it/i?) are unavoidable. In order for the scheme to yield an anomalous 
dimension of the form in Eq. (BS) these higher powers must appear along with coefficients in just the right way to 
ensure that the ln(yLt/i?) terms do not appear in 7^. This is precisely what happens for the position-scheme (jet-mass 
scheme) at two-loop order. Note that the main difference between the moment scheme and the jet-mass scheme is 
the presence of Cp terms in Jm™"™, but that there are also differences in the subleading log and constant CfPo 
terms at two-loops. In Ref. [|l9| it was proven by an explicit construction at LL order that the moment-mass has 
a consistent anomalous dimension which sums the leading logs. Thus the one-loop analysis in Ref. | p9| (which used 
the moment-mass) is fully consistent. However, beyond one-loop order it is mandatory to use the position-scheme 
definition of the jet-mass given in Eq. (57)c. 



In the position-scheme it might appear that there is a freedom in the overall normalization of Sm in Eq. (jSTpc, and in the choice of R. 
In fact to obtain a renormalon free jet-mass scheme there is no freedom in the overall normalization, there is only freedom in the choice 
of i? ||. 

^ The presence of these higher logs in -jm for the peak and moment schemes implies that these masses also do not fall into the cusp- 
anomalous dimension category, which requires an anomalous dimension of the form ■jm = R'yil'^s] + R\n{fi/R) 72[«s], i.e. with just a 
single ln{fi/R) to all orders in a^. 
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order 


peak scheme 


moment scheme 


jet-mass scheme 
e~^^ -i (5m,/ = 






C^[lnf + |] 


Cf[hif + i] 




+C7^/?o[iln^^ + |ln^ + i + f^] 
+CfCA[(|-f^)ln^ + i-4-fC3] 


C7M(4-i^)hif+8-^-2C3] 
+C^/3o[iln^# + |ln# + #-f^] 
+C^CA[(|-f^)ln^ + i-4-fC3] 




+CfC^[(i-f^)lni + ^-|^-|C3] 



TABLE I: Coefficients of the HQET counterterm 5m for different mass schemes at one and two-loop order. 



Lets extend the proof of consistency of the anomalous dimension in the position-scheme (jet-mass scheme) to the 
full non-abelian case. At the same time we will derive the very nice result that 7™ for the jet-mass scheme is entirely 
determined by the cusp-anomalous dimension. To all orders in perturbation theory, using Eq. (p7[), the jet-mass 
anomalous dimension is 
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dln/i 2 dln/i dlTi{iy) 

|) and then Eq. ( p^ ) this gives [/ij, — er^'^ /(iy)] 
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(63) 



(64) 



das{^iy)das{^l) /3[a 



da' 



where we should evaluate the final result at = R. Performing the derivatives with respect to the couplings we 
see that at any order in perturbation theory the anomalous dimension for mj{p,) is actually independent of fiy. 
Furthermore the result is given by the cusp-anomalous dimension, 7,„ = —e'^'^ RV^las^fi)]- Thus, to all orders in 
perturbation theory the jet-mass scheme, defined by c) in Eq. (57), has a consistent anomalous dimension as in 
Eq. (pq), and yields a transitive running mass, mj(/i). The final anomalous dimension equation for the jet-mass is 



dlnfi 



-e^-i?r[a,(M)], 



and is fully determined by the cusp- anomalous dimension. The all-orders solution of this equation is 
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Since is known to three-loop order we can use Eq. (E^ to obtain the running jet-mass at NNLL 
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Note that the form of the anomalous dimension in fid/dfi [mj{fi)/ R] has the same structure as that in fid/dfi [In m(/i)], 
where rn{fJ-) is the MS mass. In the remaining sections we will use the position mass-scheme and refer to it exclusively 
as the jet-mass. 



B. Relating the Jet-Mass to other Mass Schemes 

Having obtained a suitable mass definition for measurements of the top-mass from jets, we now turn to perturbatively 
connecting it to other schemes. Using the result for Smj from Table || we obtain the two-loop relation between the 
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jet-mass and pole-mass 
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This relation can be compared to other well known two-loop mass relations, such as i) between the pole-mass and 
MS-mass MM, 
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where m 
and IS-mass 
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/3o — 2/3, and as(n) here is in the (n/ -I- l)-flavor theory, and ii) between the pole-mass 
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Here Q;s(/i) is in the n^-flavor theory. Note that in Eqs. (69) and ( [70[ ) we have n/ = 5 light massless flavors, and thus 
did not write for example vacuum polarization terms depending on the b-quark mass. 

Lets imagine that the jet-mass mj{fij) is determined from a fit to massive event shapes using a scale fij T. In 
order to connect this mj{^j) to the high-energy MS-mass, we proceed as follows. First because the renormalization 
group evolution of the MS-mass does not make sense below the mass itself, we evolve the jet-mass mj{fij) up to the 
scale frit = rf%{rrit) to obtain mj{rrit) using the NNLL running result in Eq. (|67|). At this scale we then connect the 
jet and MS-masses by eliminating the pole mass from Eqs. (|6q , |69| ). Thus the two-loop relation between the jet-mass 
and MS-mass is 
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Since we take fj, = mt there is no threshold correction at the order we are working, and as{mt) in Eq. ( |7l| ) is the 
same in the n/ and (jif + l)-flavor theories. Together with Eq. ( |67| ) this formula inputs a jet-mass determined from 
production of tops far above threshold, and outputs an MS mass that can be used in other processes, such as the 
analysis of precision electroweak data. Since a high precision result for the top-mass in the IS mass-scheme can be 
determined from a threshold cross-section analysis, we also quote the two-loop conversion between the jet-mass and 
IS- mass schemes. An extra power of as is kept in the IS-scheme terms to properly ensure a renormalon free series 
At a scale fi the conversion at second order is 
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VI. RESULTS FOR THE NNLL JET FUNCTION 



In this section we present the final result for the heavy quark jet function B{s, 6m, Tt, ^), with NNLO perturbative 
corrections and a NNLL resummation of large logs. We study the numerical effect of these two-loop corrections as 
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well as of the log-resummation, including the perturbative convergence and /i-dependence of -B as a function of s, 
and in particular the stability of its peak position which is important for a top-mass measurement. At tree- level 
B{s,S'm, ^) = 5{s) and we see from Eq. ( p^ that B{s,Sm,Tt, ^) is simply a Breit-Wigner centered at s = with a 
width Tt- Beyond tree-level the jet function becomes dependent on /i and on the choice of mass-scheme through dm. 

For the cross-section (P a / dMf dM'^ in Eq. (j^) it has been proven that at any order in perturbation theory, the 
only large logs that effect the shape of the invariant mass distribution are those due to the resummation in the 
heavy-quark jet function ]l9|.^ Furthermore these large logs only exist between scales /^r T = Tt + Q^QCu/fn and 
fJ-A ^ Aqcd + rnVt/Q. The remaining large logs only modify the cross-sections normalization. The expression which 
resums all logs between the scales /xg ~ Q 3> Hm — m3>/xr — T^Ma^ Aqcd is 
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where we have defined the resummed jet function as 
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In Eqs. ( |73|j74| ) large logs are resummed by the evolution factors Uhq, Uh^, and Ub, and of these, the first two only 
affect the overall normalization. Since the scales fir and fiA differ by a factor of Q/m it is necessary to sum the large 
logs between these scales. Recall that Eq. ( [73| ) is valid for Q ^ m, which is mandatory for the top quark and antitop 
quarks to decay to well separated jets. The numerical importance of this particular resummation has already been 
demonstrated at NLL order in Ref. p9| ]. 

In the following we study the resummed jet function B{s,Sm,Tt, ^a, fJ-r) and its dependence on s and fir- In 
particular the fir dependence cancels out order-by-order in renormalization group improved perturbation theory, and 
thus the residual dependence provides a method for estimating the effect of higher order corrections to the jet 
fimction. This fir dependence cancels order-by-order between the evolutor Ub{s — s' , fi\, fir) and the fixed-order jet 
function matrix element that gives B{s' — s" , 5m, fir) in Eq. (Q). On the other hand, the dependence of the resummed 
jet function on fi\ cancels out only in the complete cross-section, where there is additional dependence on fi\ in both 
the evolution function Uh^ and the soft-function S. The analysis of the invariant mass dependence of the full NNLL 
cross-section requires constructing a consistent model for the soft-function at two-loop order, since S contains both 
perturbative and non-perturbative pieces. The procedure in Ref. can be used to carry out this analysis, but we 
leave the study of the full cross-section to a future publication. Here we focus on the resummed jet function. 

Following the strategy in appendix E of Ref. 1 19 we can obtain analytic results for the NNLL jet function even in 
the presence of the width. At NNLL order we find 
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(75) 



* In principle both the logs in the jet function and in the soft-function can modify the invariant mass distribution. However due to the 
consistency conditions discussed in Ref. ||l9[| it is always possible to exchange a summation of large logs in the soft function in favor of 
large logs in the jet function and in the hard function normalization factors. 
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The result is expressed in terms of the functions G„ = G„(s, Ff, /ia, Mr) and their s derivatives, with 
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(77) 



with H{uj) the harmonic-number function, and ip'^^\x) the fc'th derivative of the digamma function or cquivalently 
the (fc + l)'th derivative of the log of the gamma function. 

We focus our numerical analysis on two mass schemes for the jet function. In the pole scheme we take 5mi = 5m2 = 
m Eq. (H) and use a fixed pole mass m = ■mP°^'^ — 172 GeV in the formula for s in Eq. (|^). In the jet-mass scheme 
we use Smi and Sm2 from the last column of Table | with a scheme parameter R = 0.8 GeV that corresponds to a 
scale e'^'^R ~ 1.4 GeV. Here m — m,/(/ir) is the mass in the jet-scheme, and 
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The value of mj{^,Y) to be used here is obtained using the evolution equation in Eq. (|67|), running up from an input 
scale /iQ- For this scheme it is the parameter mj{po) that one will extract with future lincar-coUidcr data. In our 
analysis we take fJ,o = 2 GeV and simply fix mj{^o) = 172 GeV. We use the three-loop result for the running coupling 
everywhere. 
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where X = 1 + as{no)f3o ln(////io)/(27r) and we evolve to lower scales using the reference value as(/io = "mz) = 0.118 
with Uf — 5. Since we have systematically treated the 6-quark as massless we also ignore the 6-quark threshold in our 
coupling evolution. We also fix = 1.43 GeV, and /ia = 1 GeV. For /ip we take a central value of /ip = 5 GeV and 
consider variations about this scale in the range 3.3 GeV < fir < 7.5 GeV. Even though it may slightly underestimate 
higher-order uncertainties, we have chosen not to make the canonical choice of varying /ip up and down by a factor 
of two because of the importance of retaining the hierarchy ^r/fJ-A — Q/mj as emphasized in Ref. |^9|. 

In Fig. H we plot the resummed heavy-quark jet function B{s, dm, Ft, ^a, /^r) versus Mt. In the left panel we show 
results for the pole scheme, and in the right panel we show results for the jet-scheme. In each panel we plot tree level 
results (black dotted-line), LL results (green short-dashed lines), NLL results (blue long-dashed lines), and NNLL 
results (solid red lines). The tree results are the pure Breit-Wigner, the LL results use the tree-level B{s, Sm, fi) with 
the LL result for J7b in Eq. ([zi]), and thus correspond to just the first Gq term in Eq. (|75|). From the LL results 
we see that beyond tree-level the jet function grows a perturbative tail above the peak. The NLL results use the 
one-loop B with the NLL result for Ub and thus include the 0{as{nr)) terms in Eq. (|75|), and the NNLL curves 
use the two-loop B with the NNLL result for Ub and thus all of the terms in Eq. ([ZS]). At each of the LL, NLL, 
and NNLL orders we show three curves with /ip — 3.3, 5.0, 7.5 GeV, which are the curves from top to bottom near 
the peak respectively. Recall that in the jet scheme we fix mj{^o — 2 GeV) — 172 GeV and use the solution of the 
mass renormalization group equation in Eq. (67). Thus the conversion from Mt to s depends on the value of /ip and 
order-by-order compensates for the dependence of the residual mass terms (5mi 2(/ir)- 

Examining the LL, NLL, and NNLL results for the jet function in Fig. ^ we observe that the jet-scheme results 
in the second panel exhibit better perturbative convergence than the pole-scheme results in the first panel. This is 
true of all features, including the slope to the left of the peak, the perturbative tail to right of the peak, the peak 
location, and the peak height. Comparing the spread of the curves we see that at both NLL and NNLL order the 
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FIG. 3: The jet function, mB{s, 5m, Ft, /i) versus Mt, where s = (M^ —m'^)/m and Ft — 1.43 GeV. The left panel shows results 
in the pole-mass scheme and the right panel shows results in the jet-mass scheme. The black dotted curve is the tree-level 
Breit-Wigner, the green short-dashed curves are LL results, blue long-dashed curves are NLL, and the solid red curves are at 
NNLL order. For each of the LL, NLL, and NNLL results we show three curves with /xr = 3.3, 5.0, 7.5 GeV respectively. Other 
parameters are discussed in the text. 
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FIG. 4: Peak position Mf^^^ of the jet function versus ^y- Short-dashed results are at LL order, long-dashed are at NLL 
order, and solid are at NNLL order. Results are labeled for the pole mass-scheme (blue) and jet mass-scheme (red) . 



residual /zr dependence is smaller in the jet-scheme than in the pole-scheme. The numerical size of the residual 
scale dependence varies region by region. In the pole-mass scheme the scale dependence in the slope before the peak 
is ~ 17% at NLL and ~ 14% at NNLL, while the maximum variation near the peak is 23% at NLL and 17% at 
NNLL, and then in the tail region well above the peak it is ~ 19% at NLL and ~ 13% at NNLL. Hence, in the pole 
scheme including the NNLL results does not significantly decrease the /ir dependence. In the jet-mass scheme the 
scale dependence in the slope before the peak is ~ 6% at NLL and ~ 2% at NNLL, while the maximum variation 
near the peak is 14% at NLL and 7% at NNLL, and then in the tail above the peak it is ~ 12% at NLL and ^ 5% 
at NNLL. Thus, in the jet-mass scheme the dependence is reduced by a factor of two or more. The same level of 
improvement is observed for different values of the scheme parameter R than the value used in our analysis. 

In Fig. ^ we plot the peak position Mf'''^^ of the jet function curves, versus /ir- This figure displays the convergence 
and /^r dependence of the jet function peak position in more detail than Fig. ^. The stability of the jet function 
peak has a direct influence on the peak of the cross-section, and both are very sensitive to the value of the short- 
distance top-mass. Hence the peak-position is important to gauge the effect of perturbative corrections for the mass 
measurement. We use a wider range for /ir than that of the curves in Fig. |3[ but note that results for /ir < 3 GeV 
upset the hierarchy /ir//iA — 5 and hence can be safely ignored. In the pole-mass scheme we observe that there is 
limited sign of convergence for the peak position, although the shifts with /ir = 5 GeV at each order are still relatively 
small being ~ 230 MeV from LL to NLL order and ~ 120 MeV from NLL to NNLL order. The lack of convergence 
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of the peak-position in the pole-scheme is a reflection of the infrared renormalon in the pole-mass, which destabilizes 
perturbative predictions. We also observe from Fig. ^ that the fip dependence of the peak-position is not reduced in 
going from LL, to NLL, to NNLL in the pole-scheme. However, in the jet-mass scheme the peak location converges 
nicely from LL to NLL to NNLL, with a numerical value of Mj"'''^ = 172.099 GeV at NNLL order. Atjir = 5 GeV the 
shifts are ~ 67MeV from LL to NLL and ~ 17MeV from NLL to NNLL. Also in the jet-scheme, Fig. || shows that the 
fir dependence of the peak-position decreases from LL, to NLL, to NNLL order, with the curves becoming flatter as 
the order increases. The residual fir scale dependence of the peak position in the jet-scheme is 
where we quote the difference in M^^^^ from /ir = 3 GeV to fj,r = 10 GeV. 

Utilizing the two-loop computation with NNLL renormalization group improvement, and a jet-mass scheme with 
good renormalization group behavior, we have achieved stable results for the heavy quark jet function. However, we 
caution that the final result for the jet function is dependent on the choice of /xa, and so a more detailed phenomeno- 
logical analysis must be made only after combining the results reported here with perturbative corrections in the soft 
function to yield a ii\ independent prediction for the cross-section. In particular the size of the perturbative tail 
above the peak in the cross-section is affected by both perturbative corrections to the jet function and soft-function, 
and is strongly fi\ dependent in each of these functions individually. 



VII. CONCLUSION 



Effective field theories are an important tool for making high precision predictions for jet observables, and facilitate 
a measurement of the top mass with theoretical uncertainty less than ©(Aqcd)- Starting from the expression for 
the cross section for double top production at an e+e" collider, Eq. (|^) derived in Ref. Q, we have studied the 
properties of the heavy-quark jet function B at higher loop orders. The function B can be calculated in HQET, 



and was studied at one-loop and NLL order in Ref. |19 . Here we have performed the 2-loop computation of B to 
obtain a NNLO result. Our analysis also yielded the two- loop anomalous dimension of B, which when combined with 
the three-loop cusp anomalous dimension from Ref. Jsi] was used to obtain a renormalization group improved heavy 
quark jet function at NNLL order. Using the formulation in terms of vacuum matrix elements of Wilson lines we have 
explained precisely how the jet function is different from the heavy-quark shape function and fragmentation function 
that have been considered previously in the literature at two-loop order ^ ^ . 

The two-loop computation also allowed us to study the higher loop behavior of B and arrive at a suitable definition 
of a short distance top mass-scheme, for higher order analysis. In particular we gave a definition for a jet-mass scheme 
with nice renormalization group properties, and demonstrated that in this scheme the mass anomalous dimension 
is determined by the cusp anomalous dimension to all orders in perturbation theory. To study the properties of 
different mass-scheme definitions we exploited the fact that the heavy-quark jet function in position space fulfills the 
requirements to obey the non-abelian exponentiation theorem of Refs. p6| , |37| . In particular in the abelian limit 
{Ca, ^/ — *■ 0) the all orders result for B is simply the exponentiated one-loop result. We considered differences 
between a peak-position mass definition, a first moment mass definition, and a mass definition based on the position 
space jet function. We have checked that among these three possibilities, the peak definition and moment definition 
do not yield consistent mass renormalization group equations at LO and NLO order respectively. Thus only the 
position space definition provides a reasonable way of defining the jet-mass beyond LL order. This definition is given 
in Eq. (^)c. Relations between the jet-mass and the MS, IS, and pole masses were also given at two-loop order. The 
proof that the jet-mass is a short-distance mass, free from leading order renormalon ambiguities, is given in Ref. [ |65| . 

Our final result for the heavy-quark jet function B uses the jet-mass scheme with NNLO fixed order results and a 
NNLL resummation of large logarithms, which we refer to as the NNLL order result. The logs in this summation are 
a well defined set for a physical observable, being the only large logs that effect the shape of the top-invariant mass 
cross-section (Pa/dM^dM? . We have studied the numerical stability of B, both in terms of perturbative convergence 
from LL, to NLL, to NNLL order, and with respect to its scale dependence. In the jet-mass scheme the convergence 
of B improves by a factor of two or better in going from NLL to NNLL order. Very stable results were also obtained 
for the peak position of the heavy-quark jet function, with residual perturbative uncertainties estimated to be at the 
5 MeV level. This level of precision and stability for the jet function and its peak-position were not observed in the 
pole-mass scheme. 

Future applications of our work include the extension to complete NNLL results for the cross-section d^ a / dM^ dM^ 
by including the convolution with the soft-function and its perturbative and non-perturbative components. The use 
of position space as a convenient way of defining a top mass-scheme, could also be extended to the b-quark where 
currently a moment mass scheme, called the shape-function scheme [^, is often employed. Based on our analysis we 
expect that this shape function mass also does not have a consistent anomalous dimension beyond LL order, but that 
this can be rectified by using a modified definition of the scheme in position space. Finally, the same position space 
technique can be applied to the definition of the renormalon free gap parameter |22| , a parameter which is important 
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for combining perturbative and non-perturbative results for the soft-function for jet production. 
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APPENDIX A: TWO-LOOP GRAPHS AND RENORMALIZATION 

In this section we briefly summarize results for the two-loop jet function graphs and their renormalization factors. 
We use Feynman gauge and dimensional regularization with d = 4 — 2e. Numbering the two-loop graphs in Fig. |l| 
from left-to-right and top-to-bottom we have 

c'^f,^^ Brr^) = ^ , G. = ^ G. , (Al) 

i=l 

where we have defined G; by pulling out a common prefactor, and we let a = s + iO and t = exp(7£;)/(47r). The sum 
of terms for ^2^'° is gauge invariant. In Feynman gauge the results for the G^'s in terms of master integrals are 

Gi=ACFTFnshi, + AClh2,-CFCAh:i,, G2 - 8G|, Fo(l, 1, l)Fo(l, 1, 0) (A2) 

G3 = 2G| [Fo{l, 1, 1)] ' , G4 = -8G| F(101, 210, 100) , 

G5 = -8Gf(Gf - ]^Ca)f{\Q\, 111, 100) , Ge = -2Cf{Cf - ]^Ca)f{1{)1, 111, 101) , 

G7 = CfCa \f{111, 100, 010) - 2 F(lll, 100, 100) 4F(111, 110, 10 -1) + 2 i^(lll, 110, 000) + i^(lll, 110, 100)1 , 

Gg = CfCaF{111, 110, 100) , Gg = -4G| F(101, 110, 101) + 2CfCa -F(101, 110, Oil) , 

Gio = -CfCa [2 i^(lll, 010, 010) - i^(lll, 010, 100)] , 

Gil = {ACFTFnfh,-2CFCAf2e) [2Fo(l + e,l,l)+Fo(2 + 6,1,0)] , 

G12 = i~2CFTFnfh, + CFCAf2e) Fo{2 + e, 1, 0) , G13 = G14 = G15 = Gie = . 

Here Tp — 1/2, Ca = 3, Cf = 4/3. The Gi and the ft-^e are determined by the two-loop computation of the heavy- 
quark self-energy in Ref. js^, while the fi^ are determined from the standard sum of one- loop quark, ghost, and 
gluon vacuum polarization graphs. We have 
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The results in Eq. (A2) are given in terms of the one- loop master integral 
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(A4) 



where Nd = {in)'^^'^ fi^'^e'"^'^ , and the two-loop master integral 

F(AiA2A3,A4A5A6,A7A8A9) (A5) 

d^id^'-k Nj a(2Ai+2A2+2A3+A4+A5+A6+A7+As+A9-8) . ^jCAy+Ag+Ag) 

(277)2^ [^2]Ai[(£^_fc)2]A2[p]A3[2„.^ + a]^4[2w-(^-hfc)-fa]^=[2w-/fc-|-a]^'=[fi.^]A7[j^.£ + j^.fc]A8 [f^.fcJAg • 
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All denominator factors in square brackets in Eq. (A4) and ( |A5| ) have +iO. As written, in the light-like propagators 
[n • fc], [n ■ I + h ■ k], [fi ■ k] this prescription for -|-jO's docs not precisely match the iO's from the definition in Eq. (|44|). 
However we have checked that the difference results in scaleless integrals which are cancelled by 0-bin subtraction 
terms which are part of the definition of propagators of the coUinear fields in the jet function All 0-bin 

subtraction terms are also scaleless and therefore have not been shown explicitly. For the heavy quark jet function 
the result of adding these contributions is simply that the 1/e divergences from the integrals in Eq. (A2) are all UV. 
The general result for i^(AiA2A3, A4A5A6, AyAsAg) is not known. Several of these master integrals could be obtained 
from the calculation in Ref. |Q by shifts of variable that move a into the light-like denominators, but these are not 
the integrals needed for our analysis. We have therefore computed the master integral for the cases appearing in 
Eq. ([A^). Two cases are iterations of the one-loop master integral, -F(101, 110, 101) = — Fo(l, 1 -I- 2e, l)Fo(l, 1, 1) and 
F(101, 210, 100) = — Fo(l, 1 + 2e, l)Fp{l, 1, 0). In order to evaluate some of the remaining cases we have used the 
integration by parts technique pOl O. For simplicity we quote the results as series in e. 
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To determine the counterterms Zp'' we use the analog of Eq. (jsj) where we set Sm — and take the imaginary 
part of both sides. Since Im[mSo(s; /i)] = 5{s) this gives a simpler set of equations for the terms in the renormalized 
jet function. At one-loop we have mi?i(s, /i) = Zi{s, fi) + l'^^^ 



■2^mBi (s), and at two-loops 



mB2{s, fi) = Z2{s, fi) + i^'fi^' mBf{s) -h 2zgi fx^' mB\^'%s) + ds' Zi(s-s', /i) m^i (s') 



(A7) 



Here z„-\ = —(3q/8 enters from coupling constant renormalization. To evaluate the convolution integral term in 
Eq. (A7) we need Bi^""^ up to 0{e^), so the required ingredients from one-loop graphs are 
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where the distribution C}^ was defined in Eq. (40). From Eq. (Al) the sum of the two-loop graphs gives 
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The convolution integral required in Eq. ( A7 ) is given by 
ds' Zi{s - s', ^) t^/i2'mB5'"°(s') 
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In order to obtain this result we have used Eq. ( A8) along with Eq. (B7) from Appendix Combining the last three 
terms in Eq. (A7) the remaining l/e*^ terms must be canceled by Z2, hence uniquely fixing it. This gives 
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(All) 



Using the notation in Eq. (p5|) the counterterm consistency equations that follow from Eq. (|22|) are 
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(A12) 



Reading off the coefficients of various powers of from the results in Eqs. (A8) and (All) we can verify that they 
are all satisfied. 



APPENDIX B: RELATIONS FOR PLUS-DISTRIBUTIONS 



In the text in several places we converted between momentum space, position space, and plus-distributions arising 
from imaginary parts. Useful conversion formulas include 
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where [[p]] on the sum is the greatest integer not exceeding p. From Eq, 
Im[L"] = 6{s) , Im[ii] = , 



^.^ (2j + l)!(n-2j-l)! 



(x) hi" (x) 



a; 



(B2) 



the cases we used include 



lm[L'^] = AC-^ - Att'C^ + —5{s) , Im[i^] = -bC" + IOtt''/:^ - t:^C^ 

5 



Im[L6] = - 20Tr2£3 + eTr'^/:^ - —5{s) . 



(B3) 



Here L*-' is defined by Eq. (p^ and the distribution is defined in Eq. (Eo). The following rescaling identity is also 



useful 
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For /c > this result readily gives 



//:r/:'= = -fc£^-i-4,o<5(s) 



(B5) 



Eq. (B5) can be used to verify the expected /i-dependence at various stages. Finally, we need the convolution of two 
plus distributions. The general formula is 
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The following cases were used in our analysis 



(B6) 



w—w'—O 



1 \ (lns/M)''+^"+^ 



fc + 1 j + 1 



ds' (fi- 


s')/:«(s') 


= 2C\s) 


ds' £" (fi- 






ds' /:°(s- 







(B7) 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 



Tevatron Electroweak Working Group (2007), h ep-ex/0703034. 
M. Gruenewald, EPS-HEP, Manchester, (2007), |arXiv:0709.3744 . 
W. M. Yao at al. (Particle Data Group), J. Phys. G33, 1 (2006). 
HFAG (Heavy Flavor Averaging Group) (2007). 

I. I. Y. Bigi, M. A. Shifman, N. Uraltsev, and A. I. Vainshtein, Phys, 
A. H. Hoang, Z. Ligeti, and A. V. Manohar, Phys. Rev. D59 
A. H. Hoang, Z. Ligeti, and A. V, 



Rev. D 56, 4017 (1997), 
074017 (1999), 
Manohar, Phys. Rev. Lett. 82, 277 (1999) 



iep-ph/981123c 



bep-ph/9704245 



, , |hep-ph/9809423 

A. H. Hoang and t" Teubner, Phys. Rev. D60, 114027 (1999), [hep-ph/9904468[ 
S. W. Bosch, B. O. Lange, M. Neubert, and G. Paz, Nucl. Phys. B699, 335 (2004), |hep-ph/0402094 . 
M. Beneke, Phys. Rept. 317, 1 (1999), hep-ph/9807443. 

M. C. Smith and S. S. Willenbrock, Phys. Rev. Lett. 79, 3825 (1997), hep-ph/9612329. 

S. Fleming, A. H. Hoang, S. Mantry, and I. W. Stewart, Phys. Rev. D77, 074010 (2008), [lrep-ph/0703207 



25 



[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[32] 
[33] 
[34] 
[35] 
[36] 
[37] 
[38] 
[39] 
[40] 
[41] 
[42] 
[43] 
[44] 
[45] 
[46] 
[47] 
[48] 
[49] 
[50] 
[51] 
[52] 
[53] 
[54] 
[55] 
[56] 
[57] 
[58] 
[59] 
[60] 
[61] 
[62] 
[63] 
[64] 
[65] 
[66] 
[67] 
[68] 
[69] 
[70] 
[71] 



A. H. Hoang et al, Eur. Phys. J. direct C2, 1 (2000), hep-ph/0001286. 
A. H. Hoang, A. V. Manohar, I. W. Stewart, and T. Teubner, Phys. Rev. D65 
A. Pineda and A. Signer, Nucl. Phys. B762, 67 (2007), hep-ph/0607239. 
M. Martinez and R. Miquel, Eur. Phys. J. C27, 49 (2 003), hep-ph/02q 7315. 
A. Juste et al., ECONF C0508141, PLEN0043 (2005), |hep-ph/060lTl^ . 
A. H. Hoang and C. J. Reisser, Phys. Rev. D74, 034002 (2006), hep-ph/0604104, 
Fleming, A. H. Hoang, S. Mantry, and I. W. Stewart (2007), |0711.2079 [hep-ph] 
V. Manohar and M. B. Wise, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 10, 1 (2000) 
Neubert, Phys. Rept. 245, 259 (1994), hep-ph/9306320. 

H. Hoang and I. W. Stewart, Phys. Lett. B660, 483 (2008), |0709. 3519 [hep-ph] . 



014014 (2002), |hep-ph/0107144 . 



P. Korchemsky and S. Tafat, JHEP 10, 010 (2000), Ihep-ph/OOOTOOsT 

W. Bauer, S. Fleming, and M. E. Luke, Phys. Rev. D63, 014006 (2000), |hep-ph/0005275|. 

W. Bauer, S. Fleming, D. Pirjol, and L W. Stewart, Phys. Rev. D63, 11 4020 (2001). ^cp-ph/0011336| . 
W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D65, 05 4022 (2002'). ^iep-ph/0109045| : 
W. Bauer and L W. Stewart, Phys. Lett. B516, 134 (2001), |hep-ph/010700l| 
S. Fadin and V. A. Khoze, JETP Lett. 46, 525 (1987). 

, Beneke, A. P. Chapovsky, A. Signer, and G. Zanderighi, Phys. Rev. Lett. 93, 011602 (2 004), hcp-ph/031 2331. 
, Beneke, A. P. Chapovsky, A. Signer, and G. Zanderighi, Nucl. Phys. B686, 205 (2004), |hep-ph/0401002| . 

H. Hoang and C. J. Reisser, Phys. Rev. D71, 074022 (2005), he p-ph/0412258. 

P. Korchemsky and G. Sterman, Nucl. Phys. B437, 415 (1995), | licp-ph/9411211 . 

P. Korchemsky and G. Sterman, Nucl. Phys. B555, 335 (1999), |hcp-ph/990234l| . 

W. Bauer, A. V. Manohar, and M. B. Wise, Phys. Rev. Lett. 91, 122001 (2003), [hcp-ph/021225E|. 

W. Bauer, C. Lee, A. V. Manohar, and M. B. Wise, Phys. Rev. D70, 034014 (2004), ficp-ph/030927^ 

G. M. Gatheral, Phys. Lett. B133, 90 (1983). 

Frenkel and J. C. Taylor, Nucl. Phys. B246, 231 (1984). 

P. Korchemsky and G. Marc hesini, Nucl. Phy s. B406, 225 (1993), |hep-ph/9210281 . 
Gardi, JHEP 02, 053 (2005), |hep-ph/050i257| . 

Becher and M. Neubert, Phys. Lett. B633, 739 (2006), hep-ph/0512208. 



T 

M. Neubert (2007), arXiv:0706.2136 [hep-ph]. 
M. Neubert, Phys. Rev. D49, 3392 (1994), hep-ph/9311325. 
L L Y. Bigi, M. A. Shifman, N. G. Uraltsev, and A. L Vainshtein, Int. J. Mod 
R. L. Jaffe and L. Randall, Nucl. Phys. B412, 79 (1994), hep-ph/9306201. 

V. Manohar and I. W. Stewart, Phys. Rev. D76, 074002 (2007), [hep-ph/0605001 . 

F. Falk, M. Neubert, and M. E. Luke, Nucl. Phys. B388, 363 (1992), hep-ph/9204229. 
M. Polyakov, Phys. Lett. B82, 247 (1979). 

S. Dotsenko and S. N. Vergeles, Nucl. Phys. B169, 527 (1980). 
Brandt, F. Neri, and M. Sato, Phys. Rev. D24, 879 (1981). 

P. Korchemsky and A. V. Radyushkin, Nucl. Phys. B283, 342 (1987).^ 

Balzereit, T. Manuel, and W. Kihan, Phy s. Rev. D58, 114 029 (1998), [hep-ph/9805297 . 
Neubert, Eur. Phys. J. C 40, 165 (2005), |hep-ph/040817^ . 
J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). 
D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 
R. T. Jones, Nucl. Phys. B75, 531 (1974). 
. E. Caswell, Phys. Rev. Lett. 33, 244 (1974). 

V. Tarasov, A. A. Vladimirov, and A. Y. Zharkov, Phys. Lett. B93, 429 (1980). 
A. Larin and J. A. M. Vermaseren, Phys. Lett. B303, 334 (1993), hep-ph/9302208. 
J. Broadhurst and A. G. Grozin, Phys. Lett. B267, 105 (1991), hep-ph/9908362. 

G. Grozin (2000), hep-ph/0008300. 

Moch, J. A. M. Vermaseren, and A. Vogt, Nucl. Phys. B688, 101 (200 4), hep-ph/04031 92 



Phys. A9, 2467 (1994), |hep-ph/931235S| . 



P. Korchemsky and A. V. Radyushkin, Phys. Lett. B279, 359 (1992), |hep-ph/9203222 , 
Chay, C. Kim, Y. G. Kim, and J.-P. Lee, Phys. Rev. D71, 056001 (2005), h ep-ph/0412110. 
M. Arnesen, J. Kundu, and I. W. Stewart, Phys. Rev. D72, 114002 (2005), |hep-ph/05082il . 
Fleming, A. H. Hoang, S. Mantry, and I. Stewart, in preparation (2008). 
W. Bauer and A. V. Manohar, Phys. Rev. D70, 034024 (2004), hep-ph/0312109. 
O. Lange and M. Neubert, Phys. Rev. Lett. 91, 102001 (2003), hep-ph/0303082. 
Gray, D. J. Broadhurst, W. Grafe, and K. Schilcher, Z. Phys. C48, 673 (1990). 

Fleischer, F. Jegerlehner, O. V. Tarasov, and O. L. Veretin, Nucl. Phys. B539, 671 (1999), hep-ph/9803493. 
V. Tkachov, Phys. Lett. BlOO, 65 (1981). 

G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B192, 159 (1981). 



